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We consider the Casimir interaction between two spheres at zero and finite temperature, for both 
scalar fieids and eiectromagnetic fieids. Of particular interest is the asymptotic expansions of the 
Casimir free energy when the distance between the spheres is small. The scenario where one sphere 
is inside the other is discussed in detail. At zero temperature, we compute analytically the leading 
and the next-to-leading order terms from the functional determinant representation of the Casimir 
energy. As expected, the leading order term agrees with the proximity force approximation. The 
results for the next-to-leading order terms are new. In the limit where the radius of the outer sphere 
' goes to infinity, the results for the sphere-plane geometry are reproduced. At finite temperature, 

, the leading order term is computed and it is found to agree completely with the proximity force 

approximation in the medium and high temperature regions. For the scenario where two spheres 
^ ' are outside each other, analogous results are obtained. In the case of Dirichlet boundary conditions 

Ql on both spheres, the next-to-leading order term of the zero temperature Casimir energy is found to 

agree with that computed recently using derivative expansion. 

^ ■ PACS numbers: 12.20.Ds, 03.70.-fk., ll.lO.Wx 

^ . In recent years, the Casimir interactions between two compact objects have been under active investigation. Differ- 
^H' ent methods have been developed to compute the Casimir interactions beyond the accuracies afforded by the proximity 
force approximations. Using functional determinant or multiple scattering approach [H-Q, the Casimir interaction 
between two objects can be written in terms of the scattering matrices of each object and the translation matrices 
which relate the different bases for computing the scattering of each object. In principal, the functional represen- 
' tation of the Casimir interaction allows one to numerically compute the exact magnitude of the Casimir interaction 
' at any separation. However, the matrices involved are often infinite matrices and truncations are required, which at 
I smaller separation, pose a big challenge since the size of the matrices required for accurate computation is inversely 
[ proportional to the distance between the objects. On the analytical side, the situation is the same. It is easy to obtain 
, the large separation asymptotic expansion of the Casimir interaction from a finite number of terms in the matrices. 
However, to obtain the small separation asymptotic expansions, all the terms have to be taken into account. 

In the pioneering work a perturbation method have been developed to compute the small separation leading 
order and next-to-leading order terms of the zero terriperature Casimir interaction between a cylinder and a plane. 
This method was later extended to the sphere-plane [lol - [l^ and cylinder-cylinder [l^ configurations, as well as the 
• >— j finite temperature case Jl5J . Very recently, a new method called derivative expansion has been developed for 
] cornputing the small separation next-to-leading order term of the Casimir interaction [T6l - [l8| . While the method 
^ , of [l|, [lol - fl5j derives the asyinptotic expansion from the exact functional determinant representation of the Casimir 
■ - - ' interaction, the methods of [iGl - fTsj are completely different. In [l^, the authors considered the scalar interaction 
between an object imposed with Dirichlet boundary conditions and a Dirichlet plane. They computed the asymptotic 
expansion by performing a derivative expansion on the field in the path integral representation of the Casimir free 
energy. The leading order term coincides with the proximity force app roximation and a general expression for the 
next-to-leading order term was derived. On the other hand, in [TtI Il8j |. the next-to-leading order term is computed 
based on a postulate that the Casimir free energy have a local expansion of a certain form involving the gradient of 
the height profiles of the objects. 

In this article, we consider the Casimir interaction between two spheres, for both scalar fields and electromagnetic 
fields. This is one of the most popular configuration under consideration piitilZiiiE: 19^21]. The functional determinant 
representations of the Casimir interactions have been obtained in various works , which were used to compute 

the large separation asymptotic expansions analytically. For small separations which are more experimentally relevant. 
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it has been claimed that the first two leading terms can be computed using the derivative expansion method developed 
in ■ Since this method is based on a postulate that has not yet been proven, we find that it is necessary to compute 
the asymptotic expansions from the exact representation of the Casimir interaction. This task is undertaken in this 
paper. We will consider both the cases where one sphere is inside the other, and the two spheres are outside each 
other. 

Throughout this paper, we use units with h = c = ks = ^■ 

II. FUNCTIONAL DETERMINANT REPRESENTATION FOR THE CASIMIR FREE ENERGY 

In this section, we collect the basic formulas for the Casimir free energy between two spheres A and B with radii 
TA and tb respectively, where < rs- Let L be the distance between the centers of the spheres and let d be the 
distance between the spheres. There are two scenarios: 

• The sphere A is inside the sphere B. In this case, d = tb — rA — L. 

• The two spheres are outside each other. In this case, d = L — ta — tb- 

Starting from the representation for the zero temperature Casimir energy, the representation for the Casimir free 
energy can be obtained easily using the Matsubara formalism. The formula for the zero temperature Casimir energy 
can be derived using the multiple scattering or functional determinant approach flUol] or the mode summation approach 
[13, [l^. For a scalar field Lp, it is given by 

E^Y^' -j^f dCTrln (l - M^^iC)) ■ (1) 

Here X and Y are respectively the boundary conditions on spheres A and B. They are equal to D, N or R for Dirichlet, 
Neumann or Robin boundary conditions. In the following, when we say Robin boundary conditions, they include 
Neumann boundary condition as a special case. The trace Tr is the orbital momentum sum: 

oo oo 
r?i= — oo Z=|m| 

and is a matrix which can be written as a product of four matrices: 

^,j^XY _ rpA,y.jjABrpB.YjjBA 

The matrices U^^ and U^^ are called translation matrices. Their elements are given by 

ulpiO = i-^y^-^^^f: _(±l)''V(2^ + + 1) (o o) {L L o) 

which involve 3j-symbols. •^i"+i/2(^i) is a modified Bessel function which depends on the relative position of the 
spheres. If the sphere A is inside sphere B, 

Zv'+i/2m = ir+i/2m, 

whereas if the two spheres are outside each other, 

^i"+i/2(ei)-i^i"+i/2(ei). 

The matrices T"^^^ and T^-^, which are related to the scattering matrices, are called transition matrices. They only 
depend on the boundary conditions on the sphere A and the sphere B respectively, and they are diagonal matrices. 
For T^'-^, if the scalar field satisfies Dirichlet boundary condition (ys = on the sphere A, then the (Z, Z)-diagonal 
element is given by 
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whereas if the scalar field satisfies the Robin boundary condition dnip + aA'^ = with parameter a a on the sphere 
A, then 

' UAKi+y^{irA)+irAK[^„^{irAy 
with ua = ctA ~ 1/2- For sphere i?, if the sphere A is inside sphere i?, then 



^i + l/2(^7-B) 



respectively for Dirichlet boundary condition and Robin boundary condition with parameter as = ms + 1/2. If the 
two spheres are outside each other, then Tj^ can be obtained from the corresponding Tj^ by changing A to i?. 

For electromagnetic field, the functional determinant representation of the zero temperature Casimir energy between 
two spheres has been derived explicitly in 0, when one sphere is inside the other, and in [ol] when the two spheres 
are outside each other. It is given by 



1 r°° 

E^^/='-i^l rfeTrln(l-MW^(e)). 



(2) 



Here W and Z are respectively the boundary conditions on spheres A and B. They are equal to C or P for perfectly 
conducting or infinitely permeable boundary conditions respectively. The trace Tr is 

OO CO 

m= — oo ;=max{l,|m|} 

where the trace tr on the right hand side is a trace over 2x2 matrices: 



Here 



with 



BA 



1=0 



tttIsa} _ tAba} ^11 



A 



l{l + l)l{l + l) 



+ 1)1(1 + 1) 



The matrices Tf and are diagonal matrices given by 



t; 



TM 



If the sphere * is imposed with perfectly conducting boundary conditions, 



rp*,G,TM _ _rp*,R 



u,=l/2 



If the sphere * is imposed with infinitely permeable boundary conditions, 



T, 



*,P.TE 



T, 



u. = l/2 



t; 



.P,TM 



It is easy to see that 



pCC,T=0 _pPP,T=0 
pCP,T=0 _pPC,T=0 
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In fact, these are immediate consequences of the electromagnetic duahty. 

Using Matsubara formahsm, one can obtain the finite temperature Casimir free energy from the corresponding 
zero temperature Casimir energy by changing the integration over the imaginary frequency ^ to summation over the 
Matsubara frequencies ~ 2-K-pT . For a scalar field, this gives 

oo 

For an electromagnetic field, the matrix is replaced by the matrix M^^. In this formula, the prime ' on the 

summation over p means that the term with p = is weighted with a factor of 1/2. By taking the limit T — 0, one 
recovers the formula for the zero temperature Casimir energy. 



III. PROXIMITY FORCE APPROXIMATION 



In this section, we use proximity force approximation to find the leading asymptotic behaviors of the Casimir free 
energy and the Casimir force between two spheres when the separation between them is small. As in [l3|, we will 
consider three regions: 

(1) Low temperature region: dT ^ taT < rsT ^ 1. 

(2) Medium temperature region: dT ^ 1 ^ i'aT < rsT. 

(3) High temperature region: 1 <C <C r^T < r^T. 

First we consider the free energy density between a pair of parallel plates. Assume that the plates, plate A and plate 
are located at z = and z — d respectively, and the scalar field is imposed with the Robin boundary conditions: 

Here jSc = if the plate C is imposed with Dirichlet boundary conditions, in which case we take ac = 1- Otherwise 
Pc = 1- For this configuration, the zero temperature Casimir energy per unit area has been obtained in (23 | and 
confirmed in Q . The finite temperature Casimir free energy per unit area was computed in [25l | and it is given by 

e^id) ^:LY' r dxx\n(l- \PAX-aA){PBX-aB) 



'Cas w 27r^^^4 V {PAX + aA){liBX + aB) 

To obtain the leading asymptotic behavior of the Casimir free energy between two spheres from the proximity force 
approximation, one has to integrate {d) over one of the spheres, which we take to be sphere _B, with the distance 

d being the distance from a point on sphere B to sphere A. 

First we consider the case where sphere A is inside sphere B (see Fig. [T]). The sphere B can be parametrized by 
(0, <f): X = Tb sin0cos(/), y = tb avuO sincf), z — vb cos 9, with < 9 < tt and < 4> < 2tt. Then the proximity force 
approximation for the Casimir free energy is 

^ca^/^^ = ^l rdcl, r d9sin98lf^{h{9)) = 2TTrl f rf0 sin ^4,^ 
Jo Jo Jo 

where h{9) is the distance from the point [x, y, z) on sphere B to sphere A given by 



h{9) = ^Jrl siiY 9 + {tb cos9 - Lf -ta = + E^ - 2LrB cos9-rA. (4) 

From now on we shall assume that the spheres are eccentric, i.e., L > 0, and we are interested in the asymptotic 
behavior of the Casimir free energy when d <ti L. Making a change of variable u = h{9), we find that 



-^Cas - r' 



ri^-\-rB—rA 

■ / du(u + r^)4f^(u). (5) 

Jd 

Substituting ^ into ([S]), we have 



""^^ L f^Jd ^ 'Ji, V {pAX + aA){fiBX + aB) 
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FIG. 1: The cross section of one sphere inside the other. 



Now making a change of variables du and x x/d, we find that 



Cas 



XY,PFA_ ZbT^^, r^^^ rdxxlnU ^PAX-a^d){PBX-asd)^_,^^ 

-rA-d)j-^J^ \ 



d{rB -rA- d) Ji J^^^ V {^ax + aAd){pBX + and) 



e 



p=0 

When d <C min{rA, rs, L}, this can be approximated by 

30 fOC 



pXY.PFA ^ r-Argr r r {Pax - aAd){l3BX - asd) 



To proceed further, we need to consider the boundary conditions on the spheres. For XY = DD, /3a = = and 
we find that 

(Pax - aAd){l3BX - and) _ ^ 
{Pax + aAd){PBX + asd) 

For XY = RD, Pa ^Q,Pb = 0, and wc find that 

{Pax - aAd){PBX - aBd) _ {Pax - aAd) d->o^ Pax _ _^ 

{Pax + aAd){PBX + aBd) {Pax + aAd) Pax 

Similarly, for XY = DR, 

{Ja-i- - (Ki(l}{Pjj-i: - (isd.) d^o 



and for XY = RR, 



{Pax + aAd){PBX + aBd) ~^ 

{Pax - aAd){PBX - and) d^o^ ^ 
{Pax + aAd){PBX + aBd) 

From these, we see that in the cases of DD and RR boundary conditions, the proximity force approximation of the 
Casimir free energy gives the same leading behavior: 



EP^a^^:^!:^^' / du rfxxln(l-e-2-). 
d{rB - ta) ii J.d ^ ' 



=0 -/{pci 

The integral on the right hand side can be integrated explicitly which gives 



77PFA 

--Ad{rB-rA)\}-J^' ■ 

^ ' p=0 k=l 
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By summing over p, we have 



(7) 



where the function hs{x) is given by 

, , , 4 / coth7rA:x 1 

Differentiating with respect to d gives the proximity force approximation for the Casimir force: 



where 



\(7rfcx)2 sinh (tt/ct) (Trfcx)- 



The expressions ([7]) and ([8]) do not give us exphcitly the leading terms of the proximity force approximations to the 
Casimir free energy and Casimir force. To find the leading behaviors, we need to return to ([6]). From this expression, 
it is obvious that in the high temperature region where dT 3> 1, the dominating term is the sum of the terms with 
p = 0, which gives 



E, 



iPFA 
Cas 



tatbT J_ - 'rArBTC,R{'i) 

I,:. 



8d{rB - ta) ^ fc3 8d{rB - va) ' 



where Cni^) — X^^^i ^ ^^e Riemann zeta function. It follows that the proximity force approximation for the 
Casimir force is 

PFA rArBTCniS) 



rB -ta) 



For the low temperature and medium temperature regions where dT <^ 1, we can use the inverse Mellin transform 
formula 



1 



c-i-ioo 

dzT{z)u~ 



to transform ^ into 



M{rB - ta) -idivB - ta) J^-ioo 
Now the asymptotic behavior of when dT <C 1 can be read off from the residues of the right hand side. Using 

'z + V 



we have 

T{z)Uz)Cr{z + 3) = ^7^(^r (I) Uz)T {^-±^^ Cr{z + 3). 

Since T{z /2)(^r{z) only has poles at z = and z = 1, we find that T{z)C,r{z)(^r{z + 3) only has poles at z = 
—3, —2, —1, 0, 1. Evaluating the residues, we find that as dT <^ 1, 

pPFA T^^rArB ■K^rATBT'^ VAVBdT^ TT^rAVBd'^T'^ 

^ lUOd^rB - ta) 72{rB - ta) 2{rB - ta) ^""^ ' 90(rB -ta)' ^ ' 
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These are the only terms that are of polynomial order in dT . The remaining terms are exponentially small terms in 
dT . Differentiating with respect to d gives the proximity force approximation to the Casimir force: 

pPFA T^^TATB tatbT^ ^ ■K^rATBdT'^ 

^ 720d3(rB-r^) 2{r b - r aV''^ ' 'ih{r b - r a)' ^ ^ 

(jlOP and (jlip give the leadings terms in the medium temperature region. The first terms on the right hand sides 
are the zero temperature leading terms, and the second terms are the leading terms for the temperature corrections. 
Notice that these leading terms of the temperature corrections are finite when d — >■ 0. In the low temperature region 
where r aT < tbT <C 1, only the zero temperature leading terms are dominating. 

In the case of RD or DR boundary conditions, the proximity force approximation of the Casimir free energy gives 
the leading behavior: 



E^cl^ - , E' r du r dx X In (1 

d{rB -rA)^„ Ji Jdf, 



p=o -^1 •'d.ip 

Computing in exactly the same way as in the case of DD or RR boundary conditions, we find that 

pPFA tatbT fcC0th27rfcdr 

M{rB~rA)^, 

'^"^''^ il + Ki2dT)), 



(12) 



11520d2(rB - va) 



5760d-^[rB — ta) 



where the functions ha{x) and ga{x) are given respectively by 
720a;'* ^ f cothTrfca; 1 



ga{x) =K{x) - ^-K{x) = ^-^x^Y.{-lf-' 



coth(7rfca;) 



febi V(^^2;)^ sinh (Trfcx) (tt/cx) 



- 1. 



In the high temperature region where dT 3> 1, the leading terms for the Casimir free energy and the Casimir force 
are respectively 

pPFA irATBTCRji) 

^ 32d{rB - ta) ' 



PFA irArBT(R{3) 



Cas 



32d^rB~rA) 



When rfT < 1, 



PFA 
Cas 



Tn^rATB 



TT^rArBT^ 



TT-'rArBd^T^ 



11520d2(rs - ta) 144(rB - ta) 90(rs - r^) ' 



(14) 



PFA 7TT^rArB TT^rATBdT'^ 

^ 5760d3(rs - rA) 45(rs -ta)' ^ ' 

As before, (jl4p and (jlSp give the leading terms in the medium temperature region. The first terms on the right 
hand sides give the zero temperature leading terms, and the second terms are the leading terms of the temperature 
corrections. Notice that for the Casimir force, the leading term of the temperature correction is linear in d. The 
constant term is missing. In the low temperature region, the leading terms are the zero temperature leading terms. 
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FIG. 2: The cross section of two spheres outside each other. 



For the case where the two spheres are outside each other (see Fig. [2]) , we only integrate over the part of the sphere 
B with 

o<e <eQ^ cos-^ —. 

L 

This gives 

Ellt - 2^r| r desme£l^[h[e)), 

Jo 

where h{d) is still given by (|4]) but L = ta + + d. Making a change of variable as before, we have 

E^l^ = ^ jf^^'^^ du{u + r^)4a.(-)- 

The rest of the computations are the same as for the case where one sphere is inside the other. The results are almost 
identical to the results above, except that one has to change the factor (r^ — r^i) in the denominators to (r^ + rs)- 
In fact one finds that one would obtain the same results even if the integration is carried out over the whole sphere 
B. 

Finally, for electromagnetic field, it was well known that the Casimir free energy density for a pair of parallel plates, 
where both are perfectly conducting or both are infinitely permeable, is twice of the Casimir free energy density 
for a pair of parallel Dircichlet plates. The Casimir free energy density for a pair of parallel plates, where one is 
perfectly conducting and one is infinitely permeable, is twice of the Casimir free energy density for a pair of parallel 
plates, where one is Dirichlet and one is Neumann. Therefore, for two spheres which are both perfectly conducting or 
both infinitely permeable, the proximity force approximations for the Casimir free energy and Casimir force are twice 
of that for two spheres which are both imposed with Dirichlet boundary conditions. For two spheres where one is 
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perfectly conducting and one is infinitely permeable, the proximity force approximations for the Casimir free energy 
and Casimir force are twice of that for two spheres where one is imposed with Dirichlet boundary conditions and one 
is imposed with Robin boundary conditions. Namely, 

pPFA,CC _ pPFA,PP _ r, pPFA.DD 

pPFA.CP _ pPFA,PC _ 9 pPFA.DR 
-^Cas — ^Cas — ^-^Cas 



IV. ASYMPTOTIC BEHAVIORS OF THE ZERO TEMPERATURE CASIMIR ENERGY AND 

CASIMIR FORCE FROM EXACT FORMULAS 

In this section, we are going to compute the small distance asymptotic behavior of the zero temperature Casimir 
energy from the exact formulas given in Section |lll The method we are using is similar to that used in |l[ to find the 
asymptotic behavior between a cylinder and a plate and in (lol - [l^ to find the asymptotic behavior between a sphere 
and a plate. In [l3| . we have extended this method to the case of two cylinders. For two spheres, things are more 
involved because the M matrix is more complicated. We will explain our method in detail for the case where sphere 
A is inside sphere B. At the end of this section, we will discuss the necessary changes for the case where the two 
spheres are outside each other. 

As we mentioned in Section Hill we only consider the case where L > 0, i.e., the spheres are eccentric. Let 



rs -TA TB - TA TB - rA 

Obviously b — 1 + a. We want to find the asymptotic behavior of the Casimir energy when e <C 1. Making a change 
of variables ^ — uj/{rB — ta) and expanding the logarithm, the zero temperature Casimir energy for scalar fields ([T]) 
can be written as 



pXY,T=0 



where 



oo 

ii = \m\ 



with 



and U^^ is obtained from Uf^ by replacing li with h+i, and I'l with I'l . Here we have used the fact that 

(U U+l li \ (h h+l ii \ _ + (h+i h li \ (li+i h li\ 

[m -m OJ \0 Oj ^ ^ ' \ m J \ 00/' 

which follows from the properties of the 3j-symbols. In the following, we also have to use the fact that the 3j-symbols 

[o q) ^"""^ 1^0 o) 
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are nonzero if and only if k + U + /" and li + U+i + I'l are even. 

The small distance leading asymptotic behavior of the Casimir energy comes from the regions where Zj and u are 
large. Introduce new variables ni, . . . , n^, go, ■ ■ ■ ,qs, I'o, - ■ ■ i>o, • . • , t's and r such that 

li =1 + rii, 1 < j < s, 

h =7r('« + ^'+1) +<li = -l + ^{n-i + n-i+i) +qt, < i < s, 
la a la 

I'l - h + 2iyi = !- ~ n, + ^("-i + + + 2i^j, < i < s, ^-j^g-j 

I'l - k+i + 2i>j = - - n-i+i + ^{rii + n-j+i) + g.^ + 2;>j, < z < s, 
a la 

w = , 

ar 

where I — Iq. In terms of the new variables, the leading contribution to the Casimir energy comes from I ^ e~^, 
ni,qi,m ^ e^^/^ and VijVi^T ~ 1. The summation over h,! < i < s, is transformed to summation over Ui from 
Ui = |m| — Z to 00, which in the e ^ 1 limit can be replaced by an integration over rii from —00 to 00. The summation 
over li is transformed to summation over qi in the set 

l™l - -I - 7r-("i + + n : n = 0, 1, 2, . 

a 2a 

which in the e ^ 1 limit can be replaced by integration over qi from —00 to 00. For the summation over the 
condition |Zi — Zi| < l[' < h + h and the condition that h + U + I'l has to be even are equivalent to Vi is an integer 
satisfying max{0,Zi — Z^} < Vi < U. Therefore, in the e ^ 1 limit, the summation over I'l can be replaced by the 
summation over vi from to 00. Similarly, the summation over Z" can be replaced by the summation over Di from 
to 00. Since 

00 00 00 I 

m=~ool=\m\ 1=0 m= — I 

in the e <C 1 limit, the summations over / and m can be replaced by integrations over I and m from to 00 and from 
—00 to 00 respectively. Hence, 



(17) 



with 



In the following, we will compute the leading order term and the next-to-leading order term of the zero temperature 
Casimir energy. For this, we need to expand Mj-^y^^ up to terms of order e. First, consider the 3j-symbol 

1,1 

m —m 

Although the asymptotic expansion of this 3j-symbol has been derived in [lo| , it cannot be directly applied here since 
in our case, I'l = k — li + 2vi instead of Z" — k + k — 2vi. Nevertheless, this small problem is easy to overcome. Using 
the property of 3 j -symbols, we have 

\ ( 

m —m I \ m —m 



Similar to the case considered in flo| . we have the integral representation 

(I'l h li \ I iNm(~2i)''^''^'' [ „j /" ^ /" ^ liniB „^Ji+li~l" n ■ h+l" -L m J.\ I'-'+h-liA /I c^ 

lo m -m) ' ^i-ohm 1^9 ld(j)e cos' '6'sm- • ■(e'-0)cos' ' </) (18) 
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from the work [26|, where 



ih + my. ih - m)! ([, - m)l {k + m)l I'/l I'/l 



_{k + k + I'l + 1)! (/. + k - {h -k + i~k + k + 
Then we can proceed as in |1Q] . For the factor in front of the integral, we can use StirHng's formula 



1 



1 



In 7i! = n H — In rt — n H — In 27r 



12n 



to obtain an expansion of the form 

Oh+h+l'' - 

_2 ^I'/OUm 



exp 



m?{a + h) 



exp( Aa + A, 2 



(19) 



where here and in the following, for any X ^ A^.i and Xi^2 are respectively terms of order ^/e and e. All the expansions 
are Taylor expansions which can be done using a machine and we will not write down explicitly the expressions for 
Xi i and Xi 2- For the integral 



d6 



.e2™«cos'-+'^-'"^?sin''+'"-^ 



use the fact that sinu — u(l — v? /Q + ...), cosm = exp(— Insecu) and Insecu = + u''/12 + . . .. Introduce new 
variables p and a of order so that 



p + a 



p — aa 



one obtains an asymptotic expansion of the form 



^2ui + l roo 



da / dp a^''' 1 - 



5V 
6/ 



exp 



2im(n + a) o . t , , „ „ 
L ^ -p^ - ba^ ] {I + + B: 



i.2 



Making a change of variables p — p — ima/{bVl) and integrating with respect to p with the help of a machine, one 
obtains 



■ exp 



hi 



da a ' exp 



/ 2ima 



-ba^] 1+ C,i + C,, 



Together with (jl9p . we have the expansion 



\ V"!' — 

m -m ' ^ y 



daa'^'^' exp 



2ima 

~7r 



■ exp 



2bl 



exp^^Aa + A,2) 



(20) 



ba^][l+ C.i + C^,2 



Setting m = 0, one can integrate ([20|) with respect to a using 



da a e 



,^^„,,. ^ r(t., + i/2) ^ {2iy{)\ ^ 



These give an expansion of the form 



ii 



I'i'+h+h . 



bi Z 4 TT 4 J/ j I 



where 'Di , — Ai. 



m=0 





for = 1,2. Finally for the term 
1 



'p^2-'^--^ exp(A,i + V,^ (l + £,,i + 

4 ;a I \ /V 



(21) 



2w(l-e) 



{2k + l){2k + mi'l + l)/i;'+i/2(^(l - e)), 
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Debye asymptotic expansion of modified Bessel function gives 



where 



ry(z) =-\/l + + In 

Ul(2) = 



1 + vm^' 

1 /I 5 



VTT^VS 24(1 + ^2) 
Expanding up to terms of order e gives an expansion of the form 



V^(1-t2)i/4 \l + T, 

I T{ni + - arii + ani+i)^ rgi(nj + ra^+i - an., + an-t+i) ragf e/ \ (23) 



X exp , 

^ar 8a( 2/ 2/ ar 

X (^1 + + g^,^ exp(H,,i + ^,,2) (1 + Ji 

Here the term ^1 + Ji^ comes from the last bracket in ([2^ . with 

ar(3 - 5r2) 

— m — ■ 

which is of order e. Combining ([20)) . ([2T|) and (|23l) gives 



/ ^ T(ni + rii+i — arii + ani+i)2 Tqi{ni + Ui^i — arii -\- arii^i) Taqf el 
^ Var 8a/ 2/ 2/ ar ^ 2W 

xV^^— — / da a^"^ exp { - ba^ ] 1 + C.J + C,.2 



TV ab I 1 — T 



where 



After summing over using 



71—0 n— n— 

make a change of variable a = cr — 2(1 + r)m/(26r>/I) and integrate over a. These give 



00 



2tiI y^ + r, 

f I rijii + n^+i — arii + anj+i)^ Tqiijii + n^+i — art^ + arij+i) Taq^ el m? \ 
^ V ar 8a/ 2/ 2/ ar ^ 2Mt j 
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BA 



Interchanging and rii+i gives the corresponding expansion for t/, 

Now consider T^'^ . Debye asymptotic expansions of modified Bessel functions give 



2ui(w2) 



J • ^ _ 7]-g 



-2(/. + l/2)r,(<^2) A _ 2TOi,^^(a;2) \ 



(24) 



where 



a;6 

W2 



+ 1/2 

( ^ 1 3 7 



VTT^ V 8 ' 24(1 + z2) 
and ?7(z) and ui{z) are as before. Expanding up to terms of order e gives 

2lb aqfr qiT^Ui + n.^+i) b{ni + n^+i)'^T 



7;^'Y _(_i)a7rexp 



1 - r 



1 + r 



ar bl I Aal 

lb 1 + 



exp 



(M,1+M,2)(i+^2''), 



where y = for Y = D and y = 1 for Y = R. The term ^1 + J^^^ comes from the last bracket of the term Tf"^ in 
with 



D _ ar(3 - 5r2) _ aT{2AuB - 9 + Tr^) 

^ 12W ' 



2 126Z ' ^ 126Z 

Combining the expansions for U^^ , Uf^ and T^'^ gives an expansion of the form 



2/ r(ni + ni-j-i)"^ ar(ni — Tii+i)'^ m"^ 2£/\ /I — 



4Z 4/ bZr ar y V 1 + T 

6/ 



X ( 1 + 2Ji + J^) I dq^ exp f ) ( 1 + + 'P^,2 



(25) 



Here 

=X.,i +^,1, 



V,,2 +-^^4^4 +A^.,lM,l + +X»,2 +A^»,2 +M,2, 



M^,J ^MtJ (ui o , j = 1, 2. 

n with 

functions give 



The integration with respect to qi is straightforward. For ' , Debye asymptotic expansions of modified Bessel 



r-^i»p(2[,. + i/2l„M)(i + ^) 



iexp(2[, + l/21„(.3))(l + ?^)^) 



(26) 



where 

uja 

UJ3 



h + 1/2' 
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and ?7(z), ui(z) and uii^dz) are as before. Expanding up to terms of order e give 



(-1)^ /^l -T 

1 + r 



f - - ^ j exp(Qo + 2,^2) (1 + Js 



3 J' 



where x = if X = D, and a: = 1 if X = R. The term ^1 + comes from the last bracket of the term T^'^ in (|26p . 
with 



r(3 - 5t^ 

m 



J- 



R _ r(24MA - 9 + 7t^) 
12/ 



Combining with (j25p . we have 



2/ 



l + r 



6r 



2V7r; 



1 br 2 "^^ 2£Z\ 

"77("'- - "'+1) - ITT - — 



4/ 



blr ar J 



1 + J^Y)(l + 5,,i+5,,2), 



where 



Si,2 —'R-i,lQi,l 

a^/r 



+ 7^^.2 + Q.2, 



72, 



dgi exp 



9 ? 



J = 1,2, 



(27) 



Substitute (P7)) into PT)) . we can compute the leading order term and the next-to-leading order term of the zero 
temperature Casimir energy: 



E, 



XY,T=0 
Cas 



1 



2TrrA 



(x+y){s+l)_ 



dr- 



br 



s+l 



2^+i7r^(s + l) 7o VI - t2 7o 
m2(s + l) 2e/(s + l) 



dm j 



4/ ' '^^^ 6/r 

i=0 



s— 1 s 

1 + ^ ^ 5.,i5j,i-H^5,,2 + (s + l)j7 

i=0 j=i+l 



XY 



i=0 



We purposely separate the contribution from J'-^^ because besides the factor (— l)(^+J')(^+i), this is the only part 
that depends on the boundary conditions on the spheres. This term is independent of Ui and m. The integration over 
Hi can be performed as explained in [l], [l3| , and then the integration over m is also straightforward. After these, we 
obtain an expression of the form 



E, 



XY,T=0 
Cas 



dl I exp — 



(-^fc^) (i + r+(. + i):r-) 



(28) 



The integrations over / and r are also straightforward. For the DD case, we obtain 



E, 



DD,T=0 
Cas 



16^rf2(rs-r^)^(s + l)4 



00 -j^ 

^ (^TT) 



■K^rATB 



1 



1 + e 



1 



1 



1 



1440d2(rB - va) V TB-TA 3 
The other cases can be easily obtained from this. In the RD case, since 



3ab 
■ d d 

TA TB 



(29) 



jKD _ jDD _ jR _ jD ^ T(r2 1 + 2ua) 
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we can infer from (l28l) and (l29l) that 



E 



RD,T=0 



E 



(-1) 



s+l 



b g (-1)^+1 /•! r-Mr 

s=0 



1 

iab 
dl I exp 



2d(s + 1) 



1 



l-K^TATB 



11520d2(rB - r^) V - 3 



1 


■ d 




3 


TA 




1 


■ d 




3 







247rd ^ (s + 1)2 

s— 



40 d , 



In the DR case, 



aT(r2 - 1 + 2ub) 



Compare to the RD case, it is easy to see that 

iTT^rATB 



^ n520d^{rB - ta) 
Finally, for the RR case, we have 

RR,T=0 tatb 



tb - rA 3 



d d 
rA rB 



40 (a n 



E 



16TTd^rB-rA) ^ (s + l)4 





1 ^ 






Sab 



OO ^ 



247rd ^ (s + l) 



TT^rArB ^ d 1 1 


■ d 


d ' 


10 d 


1440^2 (rs -r^) V rB - rA 3 


rA 


rB_ 


TT^ rA 



JUA - 1) - ^ — {6UB - 1) 

TT- rs 



Next we consider the case of electromagnetic fields. From ([2]), we have 



27r(rB - r^) ^ s + 1 i, E E III E 

^ s=0 m=OZ=min{l,|m|} y=lij=min{l 



WZ,T=0 
Cas 



tr 



nM?;f.,.(-) 



where 



(30) 



Here 



/i— min{l, |m| } 



BA 



1 i 



with 



'^.('^ + 1)^.(^. + 1) 
is obtained from U^'^ by interchanging li and Z^+i 



hh 2 



T* = 



t; 



,TM 



ma;(l — e) 



k{k + l)~kik + 1) 



16 



where 



rp*,C,TM _rp*,K 



y*,P,TE _y*,R 



M, = l/2 



t: 



= 1/2 

,P,TM _ _rp*,T) 



(31) 



i" 

We need to expand A^'- up to terms of order e and Aj up to terms of order y^, which gives 



A' 



- 1- 



2lb lb 



^^/^ 



bh 



V. 



Ui is a term of order e and V is a term of order ^/e. Compare to the scalar case, it is easy to see that M^f.^^ has an 
expansion of the form: 



' exp 



T(nf-n2+i)\ (l-T 



x(-l) 
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1 + J, 




Z.TE 



1 + r 


1 + Jc 



2VttI V 4i 6/t ar 



1 + 2 J-i 1 + 5, 1 + 5, 



Z,TM 



-(1+Z^) V 
V -(1+W) 



1 + J2 





W.TE 







1 + 



W,TM 



-{l+U) V 

V -(1+w) 

(32) 



Here w = if W = C and w = 1 if W = P, z = if Z C and z = 1 if Z = P, 

'2ima 



U = 



/26t / 771^(1 + r) 

: exp ' 



Using ((3T|). we find that for fc = 2, 3 



- ba^ ]U, = - 



1 to2(1-t2) 



2blT 21252^2 



^C,TE _ 
^P,TE _ 



= 1/2 



= 1/2 
^P,TM _ 



Multiplying up the four matrices in p2p and keeping only terms up to order e in the diagonal and terms up to order 
^/e in the off-diagonal, we have 



Z,TE 







Z,TM 



1 + ^2 





W,TE 







W,TM 



-{l+U) V 



1 + J2^'™ + Js"^'™ + 2U + V^ 
~2V 



-2V 



(33) 



This term is independent of i. When we take the trace tr on Y[i=o ™f^f+i ' ^® need to take the trace of the multipli- 
cation of s -|- 1 copies of the matrix in ([33]), which up to terms of order e, is given by 



l + J,^-™ + J™ + 2z7 + V2 



tr 

^2+(s + l 



-2V 



Z,TE , ^W,TE 



-2V 

l + ji-™ + jT'^™ + 2U + V^ 



s+l 



(34) 



+ (s + 1) 



^Z,TM ^W,TM 



+ 4(s + 1)U + 2{s + 1)V'^ + 4.s(s + 1)V^ 



Substituting ([32]) and p4p into ([30]) and compare to the scalar case, we find that 



^WZ,T=0 ^ ^WZ,sc,T=0 , ^^T=0 
CcLS Ca,s Cas ' 
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where 



2'KrA 



/ ;+! ' 

2''+i7rT-fs + l 



hT 



(s + 1) h Jo 

m^{s + l) 2el{s + l)^ 



oo poo 



1=0 



blT 



dm I / drii \ 
4(s + l)z2+2(2s + l)(s + l)V2), 



CC,sc,T=0 
Cas 



pDD,T=0 pRR,T=0 



up to the next-to-leading order term, and 



^Cas — -^Cas 



UA = l/2 



E, 



DR,T=0 
Cas 



UB = l/2 



blT 



up to the next-to-leading order term. Since 

4(s + l)z7 + 2(2s + l)(s + l)V2 
it is straightforward to find that 



4{s + l) 



62 P 



Combining with the results for the scalar case, we find that 



E 



CC,T=0 
Cas 



720d2(rB - rA) 



and 



E. 



CP,T=0 



Cas 



1 



TB - TA 



TB - TA 



"1 


20' 




■ d 


d ' 


) 


3 






J A 


TB. 





80 
7^ 



rA TB 



In the following, we will briefly discuss the differences when the two spheres are outside each other. In this case, 
the small parameter is 



ta + tb 



and ^ is related to w by ^ = u)/ {rA-\-'rB)- The variables Vi and Vi are defined so that Z" — li+li — 2h'i and 1'/ — li+li — 2i>i. 
Then one has to use the integral representation given in [lo| for the 3j-symbol. The other steps are the same. We 
find that the results are similar to the case where the sphere A is inside the sphere B, one has only to change to 
—tb- This can be understood as the change in the sign of the curvature of sphere B when one change position from 
inside the sphere to outside the sphere. In the case of DD, NN and CC boundary conditions, our results agree with 
that obtained in (l7| using derivative expansion. 

As a summary, we have obtained the first two leading terms of the zero temperature Casimir energy in this section. 
For the scalar case, we have 



DD,T=0 
Cas 

RR,T=0 
RD,T=0 
DR,T=0 



TT-^rATB 



1440^2 (rs ipr^) 

n^TATB 
'~ 1440d2(rB ipr^) 

7n^rArB 
1 1520^2 (rs =Fr^) 

7n^rArB 
' 11520d^{rBTrA) 



1 ± 
1 ± 

1± 

1± 



rB T rA 
d 



1 
3 
1 

rBTrA 3 

d 1 

k - 

3 



rB T rA 
d 



1 

rB TrA 3 



d d 

— T — 

rA rB 

d d 

— T — 

rA rB 

d d 

— T — 

rA rB 

d d 

— T — 

rA rB 



20 d 20 d ,^ , 

h ^ — {SaA - 2 T ^ — {3aB - 2 , 
TT^ rA TT rB 

YTT^ TB 



(35) 
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and for the electromagnetic case, 



E 



CC,T-0 
Cas 



E. 



CP,T=0 
Cas 



720d-'irBTrA. 



1± 



1 ± 



d 



rs T rA 



"1 


20' 




■ d 




d 


3 








T 




1 


80 ■ 




' d 


T 


d 


3 ~ 











(36) 



For the terms ± and =F, the sign on the top is for the case where sphere A is inside sphere B, and the sign in the 
bottom is for the case where the two spheres are exterior to each other. As one should expect, when the two spheres 
are outside each other, all these results are symmetric with respect to rA and tb- Note that the leading terms agree 
with those obtained using proximity force approximation. 

In the scalar case, specialize the Robin conditions to Neumann conditions by setting the Robin parameters aA and 
aB equal to zero, we have 



E 



NN,T=0 



E 



ND,T=0 



E 



DN,T=0 



lMOd^{rB T rA 

Tn^rArB 
'll520d2(rB =FrA) 

Jn^rArB 
'll520d2(rB tt-a) 



1± 



1± 



1± 



rB T rA 
d 

-f 

rB T rA 

d 

h 

rB T rA 



"1 


40' 




■ d 


d ' 


3 






rA 


T — 

rB 



d d 

— T — 

rA rB 

d d 

— T — 

rA rB 



160 d 
77r2 rA 
^ 160 d 
77r2 rB 



Notice that up to the first two leading terms, we have the following relations: 



,CC,T=0 
Cas 

,CP,T=0 
Cas 



,DD,T=0 
Cas 

iND,T=0 
Cas 



,NN,T=0 
Cas ' 
,DN,T=0 
Cas 



The first relation says that the first two leading terms of the zero temperature Casimir energy between two perfectly 
conducting spheres are the sum of the first two leading terms of the zero temperature Casimir energy between two 
Dirchlet spheres and the first two leading terms of the zero temperature Casimir energy between two Neumann spheres. 
This has been claimed to be true for any two perfectly conducting bodies in [T7| . 

By taking derivative with respect to d, we find that for the zero temperature Casimir force, the first two leading 
terms are 



DD,T=0 
Cas 



■K-^rArB 



7iRR,T=0 



F 



RD,T=0 



pCC,T=0 
Cas 



^CP,T=0 
Cas 



720d3(rB T rA) 

720d3(rB T rA) 
lir^rArB 
^ 5760d3{rBTrA) 

In^rArB 
^ 5760d3{rBTrA) 
TT^rArB 
360d3(rB T rA) 
TiT^rArB 
'2880d3(rB T^-a) 



1 ± 



1 



1 ± 



2rBTrA 
1 d 
2rBTrA 
1 d 



1 ± 



1 ± 



2rBTrA 
1 d 
2rBTrA 
1 d 



1 ± 



1± 



2rBTrA 
1 d 



d d 

— T — 

rA rB 

d d 

— T — 

rA rB 

d d 

— T — 

rA rB 

d d 

— T — 

rA rB 



2rBTrA 



10 

40 " 
7^ 



TT- rA 
40 d 
77r2 rA 
40 d 
77r2 rs 

d d 

— T — 

rA rB 

d d ' 

— T — 

rA rB 



10 d 

— — (3aA 



[iUA 

(3q!b 



2 T ^— 3aB -2 



2) 

-2) 



In the case of CC boundary conditions, the next-to-leading order term has been computed numerically in [It 
reads as 



It 



F 



CC,T=0 



rArB 



Cas 



360d3(rBTM) 



1± 



2 rBTrA 



k2 d d 



where 



/ci = 1.08(±0.08), k2 = 1.38(±0.06), kz = 1.05(±0.14). 
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Our exact computation gives 

^Cas 



im<P{rB T ta) 



1 

1 ± X- 



'^rsTTA rA 



1.69 d 1.69 d 

± 



2 re 



One can see that ki agrees quite well with the exact value 1, but the errors in the fits for k2 and fcs are about 15% 
and 30% respectively. 

Finally, let us consider the case where sphere B is much larger than sphere i.e., r^i ^ r^. In the case of two 
cylinders, such scenario has been considered in [23|- Since 



1 ± 



TA 



rB T rA rs 

we find that the small separation asymptotic expansions of the zero temperature Casimir energies are given by 

TA 



E, 



DD,T=0 
Cas 

^RD,T=0 
^DR,T=0 



E 



CC,T-0 
Cas 



1440^2 
1440^2 
' 11520^2 
' 11520(^2 

TT^TA 



720^2 



1 ± 
1 ± 
1 ± 
1 ± 
1± 



TB 

rA 

TB 
TB 

^+ 

TB 

^^ 

rB 



I d ^2 d 
3 rA 3 TB 



20 
+ - 

TT 

1 80 
3^7^ 



2 (3ayi - 2 
(3aA - 2) 



± 



2 

3 TT 

2 d 



— {6aB 



d 

rB 



1 d 

3 rA 



± 



111'- 



± 

rA 3 rs 



r(3aB-2) 



E 



CP,T=o 77r3rA 



Cas 



5760^2 



1±^ 



"1 


20' 


d 


± 


"2 


20" 


3 


7r2 


rA 


3 ~ 






80 " 
7^_ 


d 
rA 


± 


2 

3-^ 


80 " 

' 7^_ 



d 

rB 



In the limit tb — > oo, we obtain the configuration of a sphere in front a plane. In this limit, (i/rs, ?'a/''b — 0. The 
results obtained above reproduce the results for a sphere in front of a plane obtained in fl^ . 



V. THE EXACT LEADING TERM OF THE FREE ENERGY 



In this section, we consider the Casimir free energy at finite temperature. We are only going to consider the 
leading term. As mentioned in Section |lll one can obtain the representation for the free energy from the zero 
temperature Casimir energy by using the Matsubara formalism, which involves changing the integration over the 
imaginary frequency ^ to summation over the Matsubara frequencies — 2'KpT . 

First consider the scalar case. From (|28|). we find that when sphere A is inside sphere B, the leading term of the 
zero temperature Casimir energy is 

Recall that 

rB - rA rAT 

Therefore, 

_ I 

Changing r back to ^ and replacing s + 1 with fc, we have 

pXY,r=o rB - (^)^ r I I 'ikd^P + [rAi? \ 
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Therefore, for the free energy, the leading term is 



Using the formula 



and integrating over we find that 

^ - ^^^^ E' E -p f-^^^^^ - ^) 



p=0 /c=l 

tatbT ^ / ir(i) 



(40) 



-iirkpdT 



tatbT _,^j,|-j._|_y)Coth27rfc(ir 



8d(rB-rA)^/ ' fc3 

Comparing to ([7]) and (If 2p . one sees immediately that this is exactly the proximity force approximation to the Casimir 
free energy. In particular, the leading terms of the Casimir free energy in the medium and high temperature regions 
agree with those predicted by the proximity force approximation. For the low temperature region, the leading term is 
the zero temperature term as predicted by the proximity force approximation. However, our method does not yield 
the low temperature leading term of the thermal correction to the Casimir free energy. Different methods have to be 
employed and we would not discuss it here. 

For electromagnetic fields, ([M)) shows that up to the leading term, 

pCC 9pDD 

These again agree with the proximity force approximations. For the case where the two spheres are outside each 
other, we also obtain the same results as those obtained by the proximity force approximations. 

If we use the method of finding the next-to-leading order term of the zero temperature Casimir energy to find the 
next-to-leading order of the Casimir free energy at the medium and the high temperature regions, we will encounter 
some difficulties resulted from the divergences in the summation over s. This signifies that either we should not Taylor- 
expand the logarithm or some of the approximations we use do not work in the medium and the high temperature 
regions. This is a complicated issue and we will not consider it here. 



VI. CONCLUSION 



In this article, we have computed the small distance asymptotic expansions of the Casimir free energy between two 
spheres from the exact representation of the Casimir free energy. We consider scalar fields with Dirchlet, Neumann or 
general Robin boundary conditions, as well as electromagnetic fields with perfectly conducting or infinitely permeable 
boundary conditions. At zero temperature, we compute the leading and the next-to-leading order terms. The results 
are summarized in (j35p and (j36p . From these, we also obtain the asymptotic expansions when one of the spheres 
becomes large. In the limiting situation where the radius of the larger sphere goes to infinity, we recover the asymptotic 
expansions for the sphere-plane geometry. At finite temperature, we obtain immediately the leading behavior of the 
Casimir free energy using the Matsubara formalism. The analytic formula we obtain agrees completely with that 
obtained using the proximity force approximation. For future work, it will be interesting to consider dielectric 
spheres. 
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